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ABSTRACT: We reconsider the possibility that Majorana masses for the three known neu- 
trinos are generated radiatively by the presence of a fourth generation and one right-handed 
neutrino with Yukawa couplings and a Majorana mass term. We find that the observed 
light neutrino mass hierarchy is not compatible with low energy universality bounds in this 
minimal scenario, but all present data can be accommodated with five generations and two 
right-handed neutrinos. Within this framework, we explore the parameter space regions 
which are currently allowed and could lead to observable effects in neutrinoless double beta 
decay, fj, — e conversion in nuclei and fi — > cy experiments. We also discuss the detection 
prospects at LHC. 
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1. Introduction 

Neutrino oscillation data [1-14] require at least two massive neutrinos with large mixing, 
providing one of the strongest evidences of physics beyond the Standard Model (SM). How- 
ever, the new physics scale responsible for neutrino masses is largely unknown. With the 
starting of the LHC, new physics scales of order TeV will become testable through di- 
rect production of new particles, so it is very interesting to explore low-energy scenarios 
for neutrino masses. Moreover, typically, these scenarios also lead to observable signa- 
tures in precision experiments, such as violations of universality, charged lepton flavour 
violating (LFV) rare decays such as £j — > ijj or fj, — e conversion in nuclei, which, being 
complementary to the LHC measurements, may help to discriminate among different mod- 
els. Regarding the fundamental question of the neutrino mass nature, Dirac or Majorana, 
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lepton-number-violating low-scale models may give additional contributions to neutrinoless 
double beta (0^2/3) decay process, shedding new light on this issue. 

On the other hand, one of the most natural extensions of the SM that has been ex- 
tensively explored in the last years is the addition of one (or more) sequential generations 
of quarks and leptons [15]. This extension is very natural and has a rich phenomenology 
both at LHC as well as in LFV processes. Moreover, new generations address some of the 
open questions in the SM and can accommodate emerging hints on new physics (see for 
instance [16] for a recent review). 

Theoretically, apart from simplicity, there are no compelling arguments in favour of 
only three families. In theories with extra dimensions one can relate the number of families 
to the topology of the compact extra dimensions or set constraints on the number of chiral 
families and allowed gauge groups by requiring anomaly cancellation. Then, one can build 
models to justify only three generations at low energies. However, one could also build other 
models in order to justify four or more generations. In the SM in four dimensions anomalies 
cancel within each generation and, therefore, the number of families is in principle free. 

From the phenomenological point of view it seems that the most striking argument 
against new generations is the measurement of the invisible Z-b oson decay width, r\ nv , 
which effectively counts the number of light degrees of freedom coupled to the Z-boson 
(lighter than which is very close to 3 [17]. However, if neutrinos from new families 

are heavy they do not contribute to Fi nv and, then, additional generations are allowed. 
Still, pairs of virtual heavy fermions from new generations contribute to the electroweak 
parameters and spoil the agreement of the SM with experiment. Global fits of models 
with additional generations to the electroweak data have been performed [18, 19] and the 
conclusion is that they favour no more that five generations with appropriate masses for the 
new particles. Although some controversy exists on the interpretation of the data (see for 
instance [20]) most of the fits make some simplifying assumptions on the mass spectrum of 
the new generations and do not consider Majorana neutrino masses for the new generations 
or the possibility of breaking dynamically the gauge symmetry via the condensation of 
the new generations' fermions; all these will give additional contributions to the oblique 
parameters and will modify the fits. Therefore, in view that soon we will see or exclude 
new generations thanks to the LHC, it is wise to approach this possibility with an open 
mind. 

From the discussion above, it seems that neutrinos from new generations are very 
different from the ones discovered up to now, since they should have a mass 10 11 times 
larger. However, this apparent difference is naturally explained within the framework that 
we are going to explore. In the SM neutrinos are massless because there are no right-handed 
neutrinos and because, with the minimal Higgs sector, lepton number is automatically 
conserved. We now know that neutrinos have masses, therefore the SM has to be modified 
to accommodate them; the simplest possibility is to add three right-handed neutrinos with 
Dirac mass terms, like for the rest of the fermions in the SM. If one then considers the 
SM with four generations (and four right-handed neutrinos), it is very difficult to justify 
why the neutrino from the fourth generation is 10 11 times heavier than the three observed 
ones. This difficulty is alleviated if right-handed neutrinos have Majorana masses at the 
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electroweak scale and the Dirac masses of the neutrinos are of the order of magnitude of 
their corresponding charged leptons [21]. Then, the see-saw mechanism is operative and 
gives neutrino masses mi ~ m^/M, 7772 ~ rrii/M, 7773 ~ m^/M, 7714 ~ m? E /M ~ uie (we 
denote by E the fourth generation charged lepton). Although with a common Majorana 
mass M at the electroweak scale it is not possible to obtain 777.3 light enough to fit the 
observed neutrino masses, this could be solved by allowing different Majorana masses for 
the different generations; but then one should explain why M2, M3 S> M4. 

Right-handed neutrinos, however, do not have gauge charges and are not needed to 
cancel anomalies, therefore their number is not linked to the number of generations. In fact, 
an extension of the SM with four generations and just one right-handed neutrino with both 
Dirac and a Majorana masses at the electroweak scale leads, at tree level, to three massless 
and two heavy Majorana neutrinos. Since lepton number is broken in the model, the three 
massless neutrinos acquire Majorana masses at two loops therefore providing a natural 
explanation for the tiny masses of the three known neutrinos [22] 1 . More generally, it has 
been shown that in the SM with ni lepton doublets, 77 # Higgs doublets and ur < right- 
handed neutrino singlets with Yukawa and Majorana mass terms there are til — tir massless 
Majorana neutrinos at tree level, of which ni — tir — max(0, n^— hh^r) states acquire mass 
by neutral Higgs exchange at one loop [24-26]. The remaining max(0,77£ — njyn^) states 
get masses at two loops. Similar extensions could be built with additional hyperchargeless 
fermion triplets, like in type III see-saw. 

In this work we reconsider the model of ref. [22], without enlarging the scalar sector of 
the SM but allowing for extra generations. The paper is organised as follows. In section || 
we summarize current neutrino data and searches for new generations. In section |3| we 
review the radiative neutrino mass generation at two loops, and show that the observed 
light neutrino mass hierarchy can not be accommodated in the minimal scenario with four 
generations. In section |3| we present a five generation example which leads to the observed 
neutrino masses and (close to tribimaximal) mixing. We introduce a simple parametrization 
of the model and explore the parameter space allowed by current neutrino data, universality, 
charged lepton flavour violating rare decays li — > £,-7 and 0^2/3 decay, as well as the regions 
that will be probed in near future experiments (MEG, fi — e conversion in nuclei). Section || 
is devoted to collider phenomenology and we summarize our results in section ||. 

2. Framework and review 

It has been well established in the last decade that neutrinos are massive, thanks to the 
results obtained with solar [1-4,11], and atmospheric [6,7,10] neutrinos, confirmed in exper- 
iments using man-made beams: neutrinos from nuclear reactors [5] and accelerators [8,12]. 

The minimum description of all neutrino data requires mixing among the three neutrino 
states with definite flavour {v e ,v tl} v T ) 1 which can be expressed as quantum superpositions 
of three massive states Vi (i=l,2,3) with masses mj. The standard parametrization of the 

1 Two-loop quantum corrections within the SM with only two massive Majorana neutrinos also lead to a 
(tiny) mass for the third one [23]. 



-3- 



leptonic mixing matrix, C/pmns> ls: 
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and Sij = sinOij. In addition to the Dirac-type phase 5, analogous to 



that of the quark sector, there are two physical phases c/>j if neutrinos are Majorana particles. 
The measurement of these parameters is by now restricted to oscillation experiments which 
are only sensitive to mass-squared splittings (Am^ = mf — m'j). Moreover, oscillations in 
vacuum cannot determine the sign of the splittings. As a consequence, an uncertainty in 
the ordering of the masses remains; the two possibilities are: 



m i < m 2 < m 3i 
7713 < mi < m-2- 



(2.2) 
(2.3) 



The first option is the so-called normal hierarchy spectrum while the second one is the 
inverted hierarchy scheme; in this form they correspond to the two possible choices of the 
sign of Arrijj = Am^ tm , which is still undetermined, while Am^ = Amj, ol is known to be 
positive. Within this minimal context, two mixing angles and two mass-squared splittings 
are relatively well determined from oscillation experiments (see table [l]), there is a slight 
hint of 613 > and nothing is known about the phases. 
Regarding the absolute neutrino mass 



scale, it is constrained by laboratory ex- 
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fects in Tritium /3-decay, which are sen- 
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The present upper limit is m Ue < 2.2 
eV at 95% confidence level (CL) [28,29], 
while a new experimental project, KA- 
TRIN [30], is underway, with an esti- 
mated sensitivity limit m Ue ~ 0.2 eV. 
However, cosmological observations pro- 
vide the tightest constraints on the ab- 
solute scale of neutrino masses, via their TaMe 1; Thg begt fit yalues q{ ^ Ught neutrinQ 
contribution to the energy density of the parameters and their la errors from [27]. 
Universe and the growth of structure. In general these bounds depend on the assumptions 
made about the expansion history as well as on the cosmological data included in the anal- 
ysis [31]. Combining CMB and large scale structure data quite robust bounds have been 
obtained: m, < 0.4 eV at 95% CL within the ACDM model [32] and J2i m i < 1-5 eV at 
95% CL when allowing for several departures from ACDM [33]. 
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Finally, if neutrinos are Majorana particles complementary information on neutrino 
masses can be obtained from 0^2/3 decay. The contribution of the known light neutrinos 
to the 0^2/3 decay amplitude is proportional to the effective Majorana mass of v e , m ee = 
| Yli m iUei\, which depends not only on the masses and mixing angles of the Upmns matrix 
but also on the phases. The present bound from the Heidelberg-Moscow group is m ee < 0.34 
eV at 90% CL [34], but future experiments can reach sensitivities of up to m ee ~ 0.01 eV [35]. 

We now briefly review the current status of searches for new sequential generations. 
Direct production of the 4 th generation quarks t' and b', assuming t' — > Wq and b' — > Wt 
has been searched in CDF, leading to the lower mass bounds m t > > 335 GeV and my > 385 
GeV [36,37]. Limits on new generation leptons, from LEP II, are weaker: mgi > 100.8 
GeV and m u i > 80.5 (90.3) GeV for pure Majorana (Dirac) particles, assuming that the 4 th 
generation leptons are unstable, i.e., their mixing with the known leptons is large enough 
so that they decay inside the detector [17]. When neutrinos have both Dirac and Majorana 
masses, their coupling to the Z boson may be reduced by the neutrino mixing angle and 
the bound on the lightest neutrino mass may be relaxed to 63 GeV [38] . While the bound 
on a charged lepton stable on collider lifetimes is still about 100 GeV, in the case of stable 
neutrinos the only limit comes from the LEP I measurement of the invisible Z width, 
m u i > 39.5 (45) GeV for pure Majorana (Dirac) particles [17]. 

Even if new generation fermions are very heavy and cannot de directly produced, they 
affect electroweak observables through radiative corrections. Recent works have shown that 
a fourth generation is consistent with electroweak precision observables [20,39,40], provided 
there is a heavy Higgs and the mass splittings of the new SU(2) doublets satisfy [40] 2 

\m t > -m b ,\ < 80 GeV, (2.5) 
\m e > -m v >\ < 140 GeV . (2.6) 

Notice, however, that a longdived fourth generation can reopen a large portion of the 
parameter space [41]. 

In addition to these phenomenological bounds one can place some upper limits by using 
perturbative unitarity, triviality and by imposing the stability of the Higgs potential at one 
loop. Typically one obtains limits of the order of the TeV [42] for degenerate lepton doublets 
and about 600 GeV for degenerate quark doublets. 

A very striking effect of new generations is the enhancement of the Higgs-gluon-gluon 
vertex which arises from a triangle diagram with all quarks running in the loop. This 
vertex is enhanced approximately by a factor 3 (5) in the presence of a heavy fourth (fifth) 
generation [39,43]. Therefore, the Higgs production cross section through gluon fusion at 
the Tevatron and the LHC is enhanced by a factor of 9 (25) in the presence of a fourth (fifth) 
generation. Thus, a combined analysis from CDF and DO for four generations has excluded 
a SM-like Higgs boson with mass between 131 GeV and 204 GeV at 95% CL [44], while 
LHC data already excludes 144 GeV < m H < 207 GeV at 95% CL [45]. From these results, 
we estimate roughly that mu > 300 GeV in the case of five generations. However, these 

2 The allowed quark mass splittings depend on the Higgs mass, according to the approximate formula 
nv -m b ,~(l + i log( n5 m g cV )) x 50 GeV from ref. [39]. 
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limits may be softened if the fourth generation neutrinos are long-lived and the branching 
ratio of the decay channel H —¥ v^v^ is significant [46]. 

Putting all together, a general analysis seems to suggest that at most only two extra 
generations are allowed [19] unless new additional physics is invoked. If extra generations 
exist, the Higgs should be heavy. Extra generation quarks should also be quite heavy and 
be almost degenerate within a generation. The constraints on new generation leptons are 
milder; charged lepton and Dirac neutrino masses should be in the range 100-1000 GeV 
and, as we will see in section 4.3 , this range will increase if neutrinos have both Dirac and 
Majorana mass terms. 



3. Four generations 

If we add one right-handed neutrino v R to the SM with three generations and we do not 
impose lepton number conservation, so that there is a Majorana mass term for the right- 
handed neutrino, a particular linear combination of u e , u^, v T , call it v' 3 , will couple to v R 
and get a Majorana mass at tree level. The other two linear combinations are massless at 
tree level but, since lepton number is broken, no symmetry protects them from acquiring 
a Majorana mass at the quantum level. In fact, they obtain a mass at two loops by the 
exchange of two W bosons (same diagram as in figure EL but with f3,f3 running in the 
loop). This leads to two extremely small neutrino masses, as desired, but there is a huge 
hierarchy between the tree-level mass, for V3, and the two-loop-level masses, for v\ and 1^2, 
therefore this possibility cannot accommodate the observed neutrino masses. 3 



w 




Figure 1: Two- loop diagram contributing to neutrino masses in the four-generation model. 

Analogously, we can extend the SM by adding a complete fourth generation and one 
right handed neutrino v R with a Majorana mass term [22,26,48,49]. We denote the new 
charged lepton E and the new neutrino v&. The relevant part of the Lagrangian is 

- 1 

£y = -£Y e e R <j) - £Y u v R 4> - -u R m R u R + H.c. , (3.1) 

3 See however [47] for a model with three generations, one right-handed neutrino singlet and two Higgs 
doublets which can accommodate neutrino masses and mixings. 
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where £ represents the left-handed lepton SU(2) doublets, eR the right-handed charged 
leptons, vr the right-handed singlet and flavour indices are omitted. In generation space 
I and cr are organized as column vectors with four components. Thus, Y e is a general, 
4x4 matrix, Y v is a general four-component column vector whose elements we denote by 
y a with a = e, fi, r, E, and niR is a Majorana mass term. The standard kinetic terms, not 
shown in eq. ([3.1|), are invariant under general unitary transformations £ — > Vj>£, eR — > V e eR 
and vr — > e ia VR. One can use those transformations, Vi and V e , to choose Y e diagonal and 
positive and also wir can be taken positive by absorbing its phase in vr. Y u is in general 
arbitrary; however, there is still a rephasing invariance in £ and eR that will allow us to 
remove all phases in Y v . 

After spontaneous symmetry breaking (SSB) the mass matrix for the neutral leptons is 
a 5 x 5 Majorana symmetric matrix which has the standard see-saw structure with only one 
right-handed neutrino Majorana mass term. Therefore, it leads to two massive Majorana 
and three massless Weyl neutrinos. From the Lagrangian it is clear that only the linear 
combination of left-handed neutrinos f 4 oc y e v e + y^v^ + y T v T + yE^E will pair up with 
ur to acquire a Dirac mass term. Thus, it is convenient to pass from the flavour basis 
(z/ e , z/^, z-v, z/g) t° a new one z/(, z/g, z/3, z/ 4 where the first three states will be massless at 
tree level and only z/ 4 will mix with vr. YiV is the orthogonal matrix that passes from 
one basis to the other we will have v a = ^2iV a i v 'i (i = a = e,fi,T,E) with 

V a 4=N a = y a / \J^2p 2/|- Since v^v^v'^ are massless, we are free to choose them in any 
combination of v e , z^, v T , ve as long as they are orthogonal to z^, i.e., ^2 a V a iN a = 
for £=1,2,3. The orthogonality of V almost fixes all its elements in terms of N a , but 
still leaves us some freedom to set three of them to zero. Following [22, 48] we choose 
V T i = Vei = Ve2 = for convenience. 

After this change of basis, we are left with a non-trivial 2x2 mass matrix for v'^ and 
vr which can easily be diagonalized and leads to two Majorana neutrinos 

1/4 = £cos6*(— 1/4 + v'l) + isin6(yR — v R ), 
1% = — sin 0(1/4 + v'l) + cos 6(vr + v R ), 



m 4) 4 = 2 m R 2 + 4m i> ^ m RJ > ( 3 - 2 ) 

where rrif) = v^J^2 i yf, with v = (<ft^), and tan 2 = m^/m^. The factor i and the 
relative signs in 1/4 are necessary to keep the mass terms positive and preserve the canonical 
Majorana condition 1/4 = z/|. If wir <C mo, we have 777.4 ~ m^, tan0 ~ 1, and we say 
we are in the pseudo-Dirac limit while when uir 3> mp, 771,4 ~ rn'o/mR and m i ~ mR, 
tan0 ~ mo/mR and we say we are in the see-saw limit. 

Since lepton number is broken by the vr Majorana mass term, there is no symmetry 
which prevents the tree-level massless neutrinos from gaining Majorana masses at higher 
order. In fact, Majorana masses for the light neutrinos, z/J,^,^, are generated at two 
loops by the diagram of figure |l|, and are given by 
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Mij = —m R m 2 D ^2 VaiV ai ml ^ VpjVptfn 2 pI a p, (3.3) 

where the sums run over the charged leptons a, (3 = e, fi, r, E while i,j = 1, 2, 3, and 

3 

I a p = J(m4,mz,m a ,mp,0) - -J(mi,m4,m a ,mp,mw), (3.4) 

with J(m4,mi,m a ,mi3,mw) the two- loop integral defined and computed in appendix |A|. 

When ran = 0, = 0, as it should, because in that case lepton number is conserved. 
Also when mo = we obtain = 0, since then the right-handed neutrino decouples 
completely and lepton number is again conserved. 

To see more clearly the structure of this mass matrix we can take, for the moment, 
the limit m e = = m T = 0; then, since we have chosen V r \ = Ve\ = Ve2 = 0, the 
only non-vanishing element in Mij is M33 and it is proportional to V^N^m^lEE- Keeping 
all the masses one can easily show that the eigenvalues of the light neutrino mass matrix 
are proportional to mf l , mf, m E which gives a huge hierarchy between neutrino masses. 



Moreover, for niE S> ^4,4 3> myy, the loop integrals in eq. (3.4) can be well approximated 
by (see appendix |A|): 

— 1 THE 

J ^"(4^| ln W (3 ' 5) 

and 

— 1 mi 

Ift/t ~ 1 tt ~ . rrz 2 In , (3.6) 

(47r) 4 2m| 7714 

leading to only two light neutrino masses, since the mass matrix in eq. ([^) has rank 2 if 
the three light charged lepton masses are neglected in I a p. The third light neutrino mass is 
generated when at least m T is taken into account in the loop integral, leading to a further 
suppression. Within the above approximation, the following ratio of V2 and 1/3 masses is 
obtained [50]: 

too 1 / to t \ / m T \ 10~ 7 
to 3 4N^\m E J \mij iVj, 

where we have taken In^mg/mi) ~ ^(ms/mj) ~ 1 and in the last step we used that 
uiE,Tng > 100 GeV. To overcome this huge hierarchy one would need very small values of 
Ne which would imply that the heavy neutrinos are not mainly ve but some combination of 
the three known neutrinos v e ,u^,u T ; but this is not possible since it would yield observable 
effects in a variety of processes, like ir — > ^lv, it — > ev, r — > euv, r — > [ivv. This requires 
that y e ^ T < 10~ 2 Ve [51,52] and then N E ~ 1. 

Therefore, although the idea is very attractive, the simplest version is unable to ac- 
commodate the observed spectrum of neutrino masses and mixings. However, notice that 
whenever a new generation and a right-handed neutrino with Majorana mass at (or below) 
the TeV scale are added to the SM, the two- loop contribution to neutrino masses is always 
present and provides an important constraint for this kind of SM extensions. 
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In the following we modify the original idea by adding one additional generation and 
one additional fermion singlet. We will see that this minimal modification is able to accom- 
modate all current data. 

4. Five generation working example 
4.1 The five generations model 

We add two generations to the SM and two right-handed neutrinos. We denote the two 
charged leptons by E and F and the two right-handed singlets by v^r and v^r. The 
Lagrangian is exactly the same we used for four generations ( |3,1[ ) but now i and e are 
organized as five-component column vectors while vr is a two-component column vector 
containing v^r and v^r. Thus, Y e is a general, 5x5 matrix, Y v is a general 5x2 matrix 
and uir is now a general symmetric 2x2 matrix. The kinetic terms are invariant under 
general unitary transformations £ — > Vpi, eR — > V e eR and v —> V v v. which can be used to 
choose Y e diagonal and positive and tur diagonal with positive elements txi^r and m^R. 
After this choice, there is still some rephasing invariance ii —> e iai £i, eiR — > e iai eiR broken 
only by Y u , which can be used to remove five phases in Y v . Therefore 

Y v =(y,i/), (4.1) 

where y and y' are five-component column vectors with components y a and y' a respectively 
(q = e, [A, T, E , F) , one of which can be taken real while the other, in general, will contain 
phases. The model, contrary to the four-generation case, has additional sources of CP 
violation in the leptonic sector. However, since at the moment we are not interested in CP 
violation, for simplicity we will take all y a and y' a real. 

Much as in the four-generation case, the linear combination v'^ oc Yl a Va^a only couples 
to v^r and the combination vL oc ^2 a y' a v a only couples to v^r. Therefore, the tree-level 
spectrum will contain three massless neutrinos (the linear combinations orthogonal to v' A and 
fg) and four heavy Majorana neutrinos. Unfortunately, since in the general case v'± and is' 5 
may not be orthogonal to each other, the diagonalization becomes much more cumbersome 
than in the four-generation case. Since we just want to provide a working example, we 
choose v'^ and u' 5 orthogonal to each other, i.e., ^2 a y a y' a = 0. This simplifies enormously 
the analysis of the model and allows us to adopt a diagonalization procedure analogous to 
the one followed in the four-generation case. 

We change from the flavour fields Vei v u,-, v ri v Ei V F to a new basis u[, u' 2 , i/' 5 where 

v[, v' 2 i ^ are massless at tree level, so we are free to choose them in any combination of the 
flavour states as long as they are orthogonal to u' A and uL. Thus, if V is the orthogonal 
matrix that passes from one basis to the other v a = £^ V a iv[ (i = 1, • • • ,5, a = e, fj,,T,E, F) 

we have V a4 = N a = y a / ^pV% K*5 = K = yUyjEpt/p, and E/3 N P N 'p = 0- The 
rest of the elements in V a i can be found by using the orthogonality of V, which gives us 12 
equations (9 orthogonality and 3 normalization conditions, because N a and iV^ are already 
normalized and orthogonal), therefore we still can choose at will three elements of V a i] for 
instance we could choose Vpi = Vf2 = Vex = 0. In this case 
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1 



m a - a = - m aR 



( Vex V e2 V e3 N e N' e \ 
V„ 2 N„ Nl 
V = V tX V t2 V t3 N t N' t . (4.2) 
V E2 V E3 N e N' e 
V V F3 N F N' F ) 

Moreover, since ^ Q N a N' a = 0, the 4x4 mass matrix of u^r, u' 5 and v^r is block-diagonal 
and can be separated in two 2x2 matrices (for v'± and v^r and u' 5 and v^r respectively) 
with the same form found in the four-generation case. Its diagonalization leads to four 
Majorana massive fields: 

V a = i COS Oai-v'a + v a) + * sin &a{v a R ~ Kr)i 

Va = - sin 9 a {y' a + v'°) + cos a (v aR + v c aR ), 

+ 4m 2 D T maR^j , (4.3) 

with a = 4, 5, tan 2 9 a = m a /m a , m 4D = VyJYla Va and m 5D = v^Y^aUa ■ 
4.2 Two-loop neutrino masses 

As in the case of four generations, the diagrams of figure [j] (now with the four massive 
neutrinos running in the loop) will generate a non-vanishing mass matrix for the three 
neutrinos u[ , v' 2 , v' 3 given by 

4 

Mij =--VE rn aR m 2 aD £ V m V aa m 2 a £ V^V^m 2 /^ , (4.4) 

m W a=4 ,5 a 

with 1^) given by ([O]) with a labeling the contribution of the 4 th and 5 th generations. 
To analyze this mass matrix first we will impose several phenomenological constraints: 

a) The model should be compatible with the observed universality of fermion couplings 

and have small rates of lepton flavour violation in the charged sector. This requires 

Ve, Vti, Vt, y' e , y'tii v'r < ve, yF,y' E , y' F - 

b) The model should fit the observed pattern of masses and mixings. A good starting 

point would be to have expressions able to reproduce the tribimaximal (TBM) mixing 
structure. 

The constraint a) together with the orthogonality condition implies that VEVE^VFyF ~ 0, 
which can be satisfied, for instance, if y F = y' F = 0, that is, v F only couples to v^r and v F 
only couples to u 5R . Then, one can define y a = y E (e e ,e^,e T ,l,0), y' a = y F (e' e ,e'^e' T ,(l : l), 
where ej and e[ are at least 0(1O~ 2 ) in order to satisfy universality constraints (see sec- 
tion 4.3.2 for more details). Thus, to order e, N a « (e e , e^, e r , 1, 0), N' a « {e! e , e' e' T , 0, 1), 



4 This pattern of couplings can easily be enforced by using a discrete symmetry which is subsequently 
broken at order e. 
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and since for i ^ 4,5 ^2 a V a iN a = ^2 a V a iN' a = 0, all the entries V a i with a = e, fj,, r, 
i = 1,2,3 can be order one. Now if we choose Vpi = Vf2 = Vei = one can see that 
Vb2,Vb3 are 0(e) while V F3 is O(e'). 

A further simplification occurs if we assume that Ve3 = 0, since in that case E only 
couples to u' 2 and and F only couples to Then, in the limit m e = = m T = the 
neutrino mass matrix My in eq. ( |4.4| ) is already diagonal and we can easily estimate the 
size of the two larger eigenvalues by neglecting the masses of the known charged leptons in 
front of rriE and mp. We find 5 



g^vn\ D vr? E the 2 g A m\D m F m F 

M 22 ~ e m m In , M 33 ~ e m 5R In . (4.5) 

"%( 47r ) 2 m i "%( 47r ) 2 m 5 



Taking 7np ~ m<5D ~ m w / 9 an d e ' ~ 10 -2 we find M33 ~2x 10 _9 m5R, therefore, to 
obtain M33 ~ 0.05 eV we need iti^r ~ 20MeV (or e' < 10~ 3 for 7715^ ~ 1 GeV). Since m§R 
and m^R control the splitting between the two heavy Majorana neutrinos we are naturally in 
the pseudo-Dirac regime unless the e's are below 10" 4 . We also see that the higher m^E^D) 
or m,E(F)> the lower the e (e') that is needed for a given m^R^Ry On the other hand, it 
is clear that the required hierarchy between M33 and M22 can be easily achieved both in 
the normal and the inverted hierarchy cases, while the degenerate case cannot be fitted 
within this scheme since the third neutrino mass is proportional to m 4 . After discussing 



the phenomenology of the model with more detail in section 4.3, we present the allowed 
regions of the parameter space in figure ||. 

Now let us turn to constraint b), that is, the light neutrino mixings. With our simpli- 
fying choices the diagonal entries of the light neutrino mass matrix are proportional to m^., 
m E , m E , whereas the off-diagonal ones are proportional to m^m E and m^m F . Therefore 
the neutrino states u[, v' 2 , u' 3 are very close to being the true mass eigenstates and the first 
3x3 elements of V, V a i, with a = e,fj,,r,i = 1,2,3 give us directly the PMNS mixing 
matrix (up to permutations). Then, by using the orthogonality conditions it is easy to 
find the structure of Yukawas that reproduce a given pattern for the PMNS matrix. Let 
us study separately the two phenomenologically viable cases, normal hierarchy (NH) and 
inverted hierarchy (IH). 

4.2.1 Normal hierarchy 

In the normal hierarchy case (mi < mi < iris), the experimental data tell us that ~ 



AmgJ ~ 0.05 eV, ~ yAm^ ~ 0.01 eV and allow for m-i <C m,2. The structures we 
have found (by choosing Vfi = Vp2 = Vei = Ve3 = 0) automatically fall in this scheme, 
since (for itie,f 3> m 4,i,5,5 ^ m w) we obtain (1711,1712,1713) oc (m^/m^,m E ,m F ). Is there 
any choice of the Yukawa couplings y a and y' a that leads naturally to some phenomeno- 
logically successful structure, for instance TBM? If we impose TBM in V a i (a = e, /j, r, 
i = 1,2,3), given the structure of N a and N' a , the orthogonality of V (at order e 2 ) imme- 
diately tells us that e e = e M = — e T = e, e' e = 0, = e' T = e' , and finally Ve2 = — ev^, 



5 Note that the position of the eigenvalues in Mij depends on the position of the zeros in Vd- The choice 



we made is very convenient to reproduce the normal hierarchy spectrum. 
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v F3 = 



—e'y/2. Therefore, a successful choice of the Yukawas will be 



y a = ys(e,e,-e, 1,0), 

y' a = y' F (0, e',e', 0,1), (4.6) 



which, keeping only terms up to order e 2 , leads to 



V 



12 1 v/3 

3 71"^"' 

1 1 ^3 







- ^ ±.±s 

V6 V3 2 v/2 V2 
J_ _L 2^ 2 J_ J_ a 





V o 



-e>/3 






-e'v/2 



— e 




1 - e' 2 / 



+ 0(e 3 



(4.7) 



Assuming that m^^ 3> 0^4 4 5 g 3> mw, we find: 



3 5 4 



m 2 = -^ml D m m m%I EE 



m 



w 



3g 4 2 2 2 m £ 

^-^ Wfi% ln-, (4.8) 



2g 

m 3 = - — 



e' 2 ml D m 5R m F I FF « f . e' 2 ml D m 5R m F In 



14 ' 



(4.9) 



and the required ratio m^jm^ ~ 5 can be easily accommodated, for instance if the fifth 
generation is heavier than the fourth one or e' > e. 



Am 2 ! I 



0.05 eV 



4.2.2 Inverted hierarchy 

In the inverted hierarchy case (7B3 < m\ < m^), we have mi ~ mi p 
and ttt-3 <C mi is allowed. Therefore now we need (mi,m2,m3) oc (m E , m' F , mf/m'g), which 
is just a cyclic permutation of the three eigenvalues. This ordering cannot be obtained 
directly with our previous choice for the zeroes in V a i, so now it is convenient to choose 
V F \ = Vf3 = Ve3 = Ve2 = instead. Following the same procedure as above, we find that 



Va = ye(-2e,e,-e, 1,0), 

y' a = y l F (e',e l ,-e l ,0,l) (4.10) 



will reproduce the desired TBM pattern, leading to 
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V3 
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l -2 C 
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+ C(e 3 ) . 



Assuming that rriE y F 3> m 4,4,5,5 ^ m VK) we get: 

3<? 4 



mi 



m 2 



2 2 z l 

e m 4D m 4 /jm £ ; In 



771E 



' (47r) 4 m^ 
2(47r) 4 m* / m 5 



(4.11) 



(4.12) 
(4.13) 

1 can be obtained by 



while the ratio of masses between the heaviest neutrinos, raxjm^ 
choosing the different parameters in their natural range. 

Thus, the model accommodates the light neutrino masses and mixings. In the next 
section we will analyse current phenomenological bounds on the mixings between the new 
generations and the first three, e, e' . 



4.3 The parameters of the model 

We have seen above that neutrino masses are proportional to e 2 m4R (or e' 2 m^) and a 
product of masses, m\ D m? E (or m^ D mp), which come from the Higgs mechanism and are 
proportional to Yukawa couplings. As discussed in section |2| the values of these masses 
cannot vary too much; perturbative unitarity requires they are smaller than about 1 TeV [42] 
while lower limits for charged leptons masses from colliders are about 100 GeV. Lower 
limits for neutral fermions are a bit less uncertain. In the case of unstable pure Dirac 
neutrinos (m a R = 0, a = 4,5) the neutrino masses are basically m a D and the lower limits 
are about 90 GeV, therefore, in that case, m a o > 90 GeV. If neutrinos have both Dirac 
and Majorana mass terms (m a R ^ 0) the masses are given by eq. ( [4.3|) and the lower limits 
are 6 m s > m a > 63 GeV, then the upper limits on m a D < ITeV automatically imply 
f^aR % 16 TeV and therefore m s < 16 TeV. More generally in figure |2| we present the 
allowed regions in the plane m a R vs m a D given the lower bound on m a > 63 GeV and the 
upper limit on m a D < 1000 GeV. We also plot the lines corresponding to m a = 200, 400, 
600 and 800 GeV. 

6 Notice that in our scenario there is a lower bound on the mixing e, in order to obtain the correct scale 
of light neutrino masses, which implies that the heavy neutrinos would have decayed inside the detector at 
LEP. 
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Figure 2: Allowed region in m a j^-m a £, given the present lower limit (m Q > 63 GeV) on the mass 
of an extra generation neutrino (a = 4, 5 refers the 4 th or 5 th generation) . 

To be definite we will take 

100 GeV < m4D,rri5£),mE,mF < 1000 GeV , 

63 GeV <m a < m aD < m- a < 16 TeV , 

with m a ma = rrv 1 ^ and m s — m a = m a R. In addition there are strong constrains from 
the electroweak oblique parameters which in the pure Dirac case require some degeneracy 
of masses, m^o ~ itie (ni^D — wif)- However, these constraints depend on the complete 
spectrum of the theory (masses of quarks and leptons from new generations and the Higgs 
boson mass) and are less certain. In fact, contributions from the splitting of masses in the 
quark sector can be compensated in part by lepton contributions with large m a R [53,54], 
which, if we do take into account the constraints set by LEP II can vary from essentially 
zero (Dirac case) to 16 TeV. 

The other parameters that enter neutrino masses are the e's, which characterize the 
mixing of light neutrinos with heavy neutrinos, and the turs, which characterize the amount 
of total lepton number breaking. The e parameters will produce violations of universality 
and flavour lepton number conservation in low energy processes. The combination of data 
from these processes will allow us to constrain both e and e'. On the other hand, to 
obtain information on the m^s we will use the light neutrino masses, which in the model 
are Major ana particles. We will also study the contributions of the heavy neutrinos to 
neutrinoless double beta decay. 

4.3.1 Lepton flavour violation processes (/i — > and \i—e conversion) 

The general expression for the branching ratio of [i — > ej produced through a virtual pair 
iy-neutrino is: 

B{^^e 1 ) = ^\8 u \\ (4.14) 
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where 

tf«, = £tf«^ (4.15) 

i 

and H is the loop function for this process [51] 

x (2x 2 + 5x - 1) 3x 3 log(x) 



H{x) 



4(x - l) 3 2(x - l) 4 



with m Xi the masses of all heavy neutrinos running in the loop and U e i and Upi their 
couplings to the electron and the muon respectively. In ( 4.14[ ) we have used the unitarity of 



the mixing matrix and neglected the light neutrino masses to rewrite the final result only 
in terms the heavy neutrino contributions. Then, as the mixings of the heavy neutrinos 
with the light leptons are different in normal and inverted hierarchy, so are the /j, — > e^f 
amplitudes generated; one can see just by inspection of the mixing matrices, eqs. (P|) and 



( 4.11 ) that in NH only the pair 1/4, 1/4 couples to both the electron and the muon, whereas 
in IH the four heavy neutrinos contribute to the process. The predicted branching ratios 
are 

NH: B{v^e 1 ) = ^Hle\ (4.16) 
IH: B^^e 7 ) = ^[H 5 e' 2 -2H,e 2 ] 2 , (4.17) 

where 

H a = cos 2 9 a H{m 2 Jm 2 v ) + sin 2 6 a H(w 2 i /m 2 v ) . 

Now since H{x) is a monotonically increasing function and > m a > 63 GeV we have 
H a > Hfaa/rriyy) > 0.09 which gives the less stringent constraint on e and e'. The 
experimental bound reads B(fi — > ej) < 1.2 x 10 -11 , and it is translated into 

NH: e < 0.03, (4.18) 

IH: |e' 2 -2e 2 |<7x 10~ 4 . (4.19) 

To see how these bounds depend on the masses of the heavy neutrinos we display in 
figure H B(/j — > ej) against the mass of the heavy neutrino 7714 in the NH case. For IH we 
expect similar results unless there are strong cancellations. We also display as horizontal 
lines present limits [17] B(fi — > e'y) < 1.2 x 10~ n and near future limits [55]. From the 
figure we can extract a conservative bound of the order of the one quoted above, e < 0.03. 

Some extra information could be extracted from r — > cy and r — > ^7. Thus, from 
B(t ej) < 3.3 x 10~ 8 we obtain e < 0.3 in the case of NH and \e' 2 — 2 e 2 | < 0.08, limits 
that show exactly the same dependence on e and e' as the one obtained in fi — > ej, but 
which are roughly one order of magnitude worse. From B(t — > ^) < 4.4 x 10~ 8 , although 
the bounds are of the same order of magnitude as for B(t — > ej), we obtain different 
combinations of e's, |e' 2 - e 2 | < 0.09 for NH and e' 2 + e 2 < 0.09 for IH. 

Another very interesting process which gives information on e is fi-e conversion in 
nuclei. From present data [17] one obtains bounds similar to the limit obtained from 
[i — > e'y. However, there are plans to improve the sensitivity in fi-e conversion in 4 and 
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Figure 3: B(fi — > ey) against 7714 for different values of e in the NH case. We also display present 
and future limits on B(fi ej) as horizontal lines. 



even 6 orders of magnitude [56], therefore we expect much stronger bounds in the future 
coming from fi—e conversion. Strong correlations between both processes exist, as can be 
seen in [51]. 

4.3.2 Universality bounds 

New heavy generations that couple to the observed fermions can potentially lead to viola- 
tions of universality in charged currents because of the "effective" lack of unitarity in the 
mixings when the heavy generations cannot be produced. Data from neutrino oscillation 
experiments can also be used to constrain deviations from unitarity of some of the elements 
of the leptonic mixing matrix [57], however, in our scenario they lead to weaker bounds 
than the ones obtained here. 

There are different types of universality bounds which constrain the mixings of light 
fermions with new generations: 

• Lepton-hadron universality. One compares weak couplings of quarks and leptons 
using muon decay and nuclear f3 decay, which are very well tested. In our case this 
involves mixings both in the quark and lepton sectors and they are not useful to test 
individually the lepton mixings we are interested in. 

• Relations between muon decay, mz , rn\v an d the weak mixing angle sin 2 9yy . These 
are very well-determined relations in the SM, and in our case they are modified because 
the heavy neutrinos cannot be produced in ordinary muon decay. Unfortunately, these 
relations depend strongly on the p parameter, which receives contributions from the 
Higgs and very large contributions from the heavy fermions of the new generations. 
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Therefore, although these type of relations could be used to set bounds on the e's, 
they would depend on other unknown parameters. 



• Ratios of decay widths of similar processes. The bounds obtained from this type of 
processes are very robust because most of the uncertainties cancel in the ratios. We 
will only consider the most precise among these ratios, which are well measured and 
can be computed accurately [58-60]: 

— r(vr — > ev) 
„ T(t — > evv) 

Rr^e/r^ = p ~ - , 4.21 

__ T(r — > evv) 

Rr^e/a^e = TT, ; - \ = 5 T-}e — , (4.22) 

— > evv) T T 

R ^l^ = ffo -+ evv) ~ Br ^- T > (4 ' 23) 

where B T ^f = T(t — > fvv)/T(r — > all) is the branching ratio of the tau decay to the 
fermion /, and r/ = 1/T(/ — > all) its lifetime. In our model there are corrections to these 
ratios because v e , v^ and v T have a small part of v A 4 and v 5 5, which are heavy and cannot 
be produced. This leads to an additional violation of universality which depends on the 
mixings of v e , v^ and v T with v 4 ^ and f 5j g. For R^^/^^.^, and using the V a i in (JOj) and 



(|4.1l[) , we find that 

Rn^e/n^ = \ Vel \ 2 + | Ve2 1 2 + I Ve3 f = { 1 + ^ NH 

*W,r-„ l^i| 2 + |^ 2 | 2 + |^ 3 | 2 1 1 " 3e 2 IH 



(4.24) 



R T ^e/r^^/ Rf^ie/r^p tests exactly the same couplings, therefore the result is the same 
as in (fO|). 



Rr^e/fi^e gives a different information because it tests r/fi universality; however, we 
find that for both NH and IH R T ^ e/ 'R^ e/ = 1, and since it is independent of the 
e's, this process does not give any further information. This is a consequence of our choice 
for the Yukawa couplings 7 , which, up to signs, are equal for the r and [i neutrinos. 

Finally for R T _^^i^ e , using our mixing matrices, we find 

Rr^/^e |Kl| 2 + |K 2 | 2 + |^r 3 | 2 J 1 - ^ NH 

\Ve,\ 2 + \Ve2\ 2 + We,\ 2 ^ 1 + 3e 2 IH [ ^ 

which gives exactly the inverse combinations of those obtained from e/fi universality tests. 



7 Which, in turn, is a consequence of the TBM structure we wanted to reproduce. 



Observable 


12 [17] 


R SM 


R/R SM 




(1.230 ±0.004) x 10" 4 


(1.2352 ±0.0001) x 10" 4 


0.996 ± 0.003. 


Rr— >e/r— >/i 


1.028 ± 0.004 


1.02821 ±0.00001 


1.000 ± 0.004 


T — 11— >e 


(1.31 ±0.06) x 10 6 


(1.3086 ±0.0006) x 10 6 


1.001 ±0.004 



Table 2: Relevant universality tests 



Therefore, if we combine the three results, R^^/^JR^^^, R T ^ e / T ^/ 'R™ e/T ^^ 
and {R T ^n/ Li^re/ R^^^/ ^e)^ 1 an( i use the data collected in table (g), we obtain 



1 + e' 2 NH 

0.998 ± 0.002 = { 1 _ 3e 2 m . ( 4 - 26 ) 



which translates into the following upper 90% C.L. limits on e' and e 



NH : e' < 0.04, (4.27) 
IH : e < 0.04. (4.28) 

Notice that although in the IH case we have more sensitivity than in the NH case because 
of the factor of 3 we finally obtain similar limits in the two cases. This is because in the 
NH case the deviation from 1 obtained in the model is always positive while the present 
measured value is slightly smaller than 1 (in both cases we used the Feldman & Cousins 
prescription [61] to set 90% C.L. limits). 

Now we can use all data from LFV and universality and conclude that in the NH case 
we have e < 0.03, basically from \i — > ej, and e' < 0.04, basically from universality tests. 
In the IH case we obtain that, except in a narrow band around e' 2 ~ 2e 2 , e < 0.02 and 
e' < 0.03 basically from \i — > if e' 2 ~ 2e 2 there is a cancellation in \i — > e^f but still one 
can combine these data with the universality limits to obtain e < 0.04 and e' < 0.06. 

4.3.3 Neutrinoless double beta decay (0v2(3) 

As commented above mm^m^R are the relevant parameters which encapsulate the non- 
conservation of total lepton number and they control, together with the e's, the neutrino 
masses. Therefore, it would be useful to have additional independent information on these 
parameters. The most promising experiments to test the non-conservation of total lep- 
ton number are neutrinoless double beta decay experiments. The standard contribution, 
produced by light neutrinos, to 0^2/3 has largely been studied (for a recent review see for 
instance [35]) and, given the expected future sensitivity, m ee = 0.01 eV, it will be very 
difficult to see it unless the neutrino spectrum is inverted or degenerate. However, if heavy 
neutrinos from new families are Majorana particles, they lead to tree-level effects in neu- 
trinoless double beta decay [62], while, in our scenario, light neutrino masses are generated 
at two loops; thus, in principle, it is possible that these new contributions dominate over 
the standard ones. 
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The contribution of new generation heavy neutrinos (with mass larger than about the 
proton mass, m p ) to the rate of neutrinoless double beta decay can be written in terms of 
an effective mass, 

(Mn)- 1 = J2ulM-\ (4.29) 

a 

where U ea is the coupling of the electron to the left-handed component of the heavy neutrino 
a. The non-observation of neutrinoless double beta decay implies that [63] 

(M N ) > 10 8 GeV . (4.30) 

In the case of NH the electron only couples to z/4 and v~i (see eqs. ( |4.7D and ( [4. 3D ), thus 

, „ cos 2 9 sin 2 n itiar 

(Mm)- 1 = 2e 2 ( ) = e 2 ^ , (4.31) 

vtla m 4D 



and using ( 4.30| ), we get: 

m iR e 2 /ml D < lO^GeV" 1 . (4.32) 

This is the same combination of the relevant parameters (m^ue 2 ) that appears in ([4.8|) for 
rri2, therefore we can use neutrino data to set bounds on the heavy neutrino contribution 
to neutrinoless double beta decay written in terms of the effective mass. We obtain 

( M n) = 9 (A ,4 4™* £ 2 x 10 11 GeV , (4.33) 

2TO 2 (47r) 4 77% 

where we have used m,2 ~ 0.01 eV, typical values for TJI40 ~ txie ~ 100 GeV and h^m^/mj) r 



1. This is far from present, eq. 4.30, and future sensitivities. 



In the case of IH, the effective mass is given by 

(Mn)- 1 = 4e 2 m 4R /ml D + e' 2 m 5R /m 2 5D , (4.34) 

leading also to unobservable effects in 0v2{3 decay. 

To summarize all phenomenological constraints on the model, in figure |] we show in blue 
the allowed region in the e — mm plane, which leads to M33 ~ 0.05 eV varying the charged 
lepton masses itie {tuf) and the Dirac neutrino masses m^o ("I5Z)) between 100 GeV-1 TeV, 
and imposing the LEP bound on the physical neutrino mass, 7714 > 63 GeV. We also plot 
the present bounds on the mixings e (e') from fj, — > &y and future limits from /i-e conversion 
if expectations are attained. 



5. Collider signatures 

As we mentioned before, the LHC offers a unique opportunity to discover (or exclude) 
new sequential generations of quarks and leptons. For instance, with 1 fb _1 at 7 TeV, the 
exclusion bound on b' would reach 500 GeV via b' — > Wt decay channel, close to the partial 
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Figure 4: Parameter space that predicts the right scale for heavy and light neutrinos (blue region 
between the curves). As a comparison we also present the current bound from /x — > ej and future 
limits from fi—e conversion experiments. 

wave unitarity bound. Even if the t' and b' are too heavy to be seen directly, their effects 
may be manifest at LHC, since they induce a large gg — > ZZ signal [64]. See also [41] for 
prospects of detecting very long-lived fourth generation quarks, i.e., in the case of extremely 
small mixings with the lighter three generations. 

Regarding the lepton sector, the standard searches for a fourth generation have to 
be restricted to the parameter space which leads to the correct light neutrino mass scale, 
depicted in figure |||. The expected signatures depend on the nature (Dirac or Majorana) of 
the neutrinos, which are generally assumed to be the lightest states. 

Most theoretical analysis of fourth generation Majorana neutrino at hadron colliders 
have focused on the process qq' — > —¥ v^i^, where the fourth generation neutrino is 
produced in association with a light charged lepton [65-68]. Subsequently, if neutrinos are 
Majorana, they will decay through 1/4 — > W^l^, leading to the low-background like-sign di- 
lepton signature in half the events. However in our model the cross section for this process 
is suppressed both by the mixing of the extra generations with the first three and by the 
small Majorana masses ui^r^r, much as in the neutrinoless double beta decay discussed 
above, so it will not be observable at LHC for the parameter range that reproduce the 
correct scale of light neutrino masses. 

Alternatively, the lighter neutrinos can be pair-produced via an s-channel Z boson, 
qq — > Z — > ujuj (I, J = 4,4) [69]. Although the W production has a higher cross section 
than the Z at hadron colliders, and the mass reach is enhanced when only one heavy particle 
is produced, if the mixing angle between the extra generations and the light ones is less 
than about 10 -6 the neutrino production rates in the W channel are so suppressed that 
they are unobservable [65]. However the rate of heavy neutrino pair- production via a Z 
boson is independent of this mixing, becoming the dominant production mechanism in the 
small mixing regime. Moreover, if the mass difference between 1/4 and 1/4 is at least 1 GeV, 
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and the mixing so small that the decay v\ — > v±Z always dominates, the above processes 
also lead to like-sign leptons in half of the events 8 . See ref. [69] for a detailed study of the 
Tevatron dataset potential to exclude (or discover) fourth generation neutrinos with both, 
Dirac and Majorana masses, up to 150-175 GeV, depending on the mixing. For the LHC, 
only the pure Majorana case has been studied in ref. [70]. According to them, the LHC at 
y/s = 10 TeV with 5 fb _1 could expect to set a 95% CL mass lower limit of m^r > 300 GeV 
or report 3<7 evidence for the vt± if m Ui < 225 GeV. We expect a similar sensitivity in 
our model, in the region 7714 — > 1 GeV and small mixing (e,e' < 10~ 4 ) i.e., somehow 
complementary to the one probed in LFV processes. See also [71] for an evaluation of the 
LHC discovery potential for both Majorana and Dirac type fourth family neutrinos in the 
process pp — > Z/H — )• — > WfiWfi. 

Searches for fourth generation charged leptons at the LHC have been studied in [72], 
also in a general framework with Dirac and Majorana neutrino masses, and assuming that 
the neutrino 1^4 is the lightest fourth generation lepton. For charged leptons with masses un- 
der about 400 GeV, the dominant production channel is charged lepton - neutrino, through 
the process qq' — > — > v^E . The neutrino v± can only decay to U4 — > W£, and being 
Majorana it can decay equally to W~£ + and W + £~ . Therefore when a pair of fourth gen- 
eration leptons are produced, we expect the decay products to contain like-sign di-leptons 
half of the time. The sensitivity study for this process in events with two like-sign charged 
leptons and at least two associated jets shows that with y/s = 7 TeV and 1 fb _1 of data, 
the LHC can exclude fourth generation charged lepton masses up to 250 GeV. It would 
be interesting to study the parameter space in our model that would lead to this type of 
signals. 

In the above searches, it was assumed that the lightest neutrino decays promptly. How- 
ever, if the mixing of the lightest fourth generation neutrino with the first three generation 
leptons is e < 10 -7 its proper lifetime will be T4 > 10~ 10 s. The decay length at the LHC 
is given by d = ficyr^ ~ 3cm(r4/10 _10 s)/37, thus for T4 > 10~ 10 s the fourth neutrino will 
either show displaced vertices in its decay or decay outside the detector, if d > O(m), which 
is a typical detector size. In our scenario, such a tiny mixing is only compatible with large 
Majorana masses, ~ 1 TeV (see figure f|), far from the pseudo-Dirac case. Searches 
for Majorana neutrinos stable on collider times have been discussed in [73], where it is 
proposed to use a quadri-lepton signal that follows from the pair production and decay of 
heavy neutrinos pp — > Z — > v^v^ — > v^v^ZZ , when both Z's decay leptonically. The final 
state is thus 4£ plus missing energy. For 30 fb _1 of LHC data at 13 TeV, 1^4 masses can be 
tested in the range 100 to 180 GeV, and 1^4 masses from 150 to 250 GeV. 

Finally, if the lightest fourth generation lepton is the charged one, there is a striking 
signal which to our knowledge has not been studied in the literature: lepton number vio- 
lating like-sign fourth generation lepton pair-production, through qq' — > — > E < ^44 — > 
or via W fusion, qq W ± W ± q'q' -)• E ± E ± jj These processes are not sup- 
pressed by the small mixing with the first three generations, so in principle they could 

8 In the exact Dirac limit, v± must decay to W£ and the different contributions to same sign di-lepton 
production cancel, since the Dirac neutrino conserves lepton number. However, as far as always decays 
to 1/4Z there is no interference amplitude, and same-sign di-lepton decays are unsuppressed. 
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be observable in our scenario. Depending on the charged lepton lifetime, they will decay 
promptly to same-sign light di-leptons, show displaced vertices or leave an anomalous track 
of large ionization and/or low velocity. A detailed phenomenological study would be very 
interesting, but it is beyond the scope of this work. 

6. Summary and conclusions 

We have analysed a simple extension of the SM in which light neutrino masses are linked to 
the presence of n extra generations with both left- and right-handed neutrinos. The Yukawa 
neutrino matrices are rank n, so if we do not impose lepton number conservation and allow 
for right-handed neutrino Majorana masses, at tree level there are 2n massive Majorana 
neutrinos and three massless ones. In order to obtain heavy neutrino masses above the 
experimental limits from direct searches at LEP, the Dirac neutrino masses should be at the 
electroweak scale, similar to those of their charged lepton partners, and the right-handed 
neutrino Majorana masses can not be too high (of order 10 TeV at most). The three 
remaining neutrinos get Majorana masses at two loops, therefore this framework provides 
a natural explanation for the tiny masses of the known SM neutrinos. On the other hand, 
it should be kept in mind that the two-loop contribution to the neutrino mass matrix is 
always present in this type of SM extensions, therefore the experimental upper limit on the 
absolute light neutrino mass scale leads to a relevant constraint which has to be taken into 
account. 

We have shown that the minimal extension with a fourth generation can not fit si- 
multaneously the ratio of the solar and atmospheric neutrino mass scales, Am^/Am^, 
the lower bound on the heavy neutrino mass from LEP and the limits on the mixing be- 
tween the fourth generation and the first three from low energy universality tests. Then, 
there are two possibilities: either enlarge the Higgs sector [47] or consider a five generation 
extension. In this work we have analyzed the second one, while the first will be studied 
elsewhere [74] . Notice that five generations are still allowed by the combination of collider 
searches for its direct production, indirect effects in Higgs boson production at Tevatron 
and LHC, and precision electroweak observables [19], provided the Higgs mass is roughly 
mil > 300 GeV. However they will be either discovered or fully excluded at LHC, making 
our proposal falsifiable in the very near future. 

Given the large number of free parameters in a five generation framework (10 neutrino 
Yukawa couplings, 2 charged lepton masses and 2 right-handed neutrino Majorana masses), 
we have considered a very simple working example assuming that i) the linear combinations 
of left-handed neutrinos that get Dirac masses at tree level are orthogonal to each other, 
ii) each extra generation left-handed neutrino couples only to one of the two right-handed 
SM singlet states and iii) each extra generation charged lepton couples only to one linear 
combination of the (tree-level) massless neutrinos. Then, we are left with 2 neutrino Yukawa 
couplings yE,UF, 2 charged lepton masses mp, mp, 2 right-handed neutrino Majorana 
masses m a R,a = 4, 5, which characterize the amount of total lepton number breaking, and 
two small parameters e, e' which determine the mixing among the first three generations 
and the new ones. Moreover, at leading order the two-loop neutrino masses 777-2, ™-3 depend 
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only on e, 7TI4R, yg, Trig and e', 7715^, yjr, tuf, respectively (see eq. (|4.5|) ). Even in this over- 
simplified case we are able to accommodate all current data, including the observed pattern 
of neutrino masses and mixings both for normal and inverted hierarchy spectrum. A definite 
prediction of the model (independent of the above simplifying assumptions) is that the three 
light neutrinos can not be degenerate. 

We have explored the parameter space regions able to generate the correct scale of 
neutrino masses, ~ 0.05 eV. We find that for typical values of rriE,m4D (ttif, 771515) at the 
electroweak scale, we need e 2 m4R {(-^m^R) < 1 keV to obtain the atmospheric mass scale 
(see figure @j). 

We have also studied the current bounds on the mixing parameters e and e' from the 
non-observation of LFV rare decays £ a —> £^7, as well as from universality tests. All of 
them are independent of the Majorana masses mm, since they conserve total lepton number. 
Depending on the light neutrino mass spectrum (normal or inverted), the strongest bounds 
come from /i — > e'j and from universality tests in n decays. Combining the information from 
both processes we can set independent limits on e and e' which being quite conservative are 
of the order of the few percent, e < 0.03 and e' < 0.04. 

Finally, we have analysed the phenomenological prospects of the model. With respect 
to LFV signals, future MEG data will improve the limits on the e's by a factor of about 3 
while, if expectations from fi-e conversion are attained the limits on the e's will be pushed to 
10 -3 . This region of observable LFV effects corresponds to the pseudo-Dirac limit, ra a n < 
1 GeV, i.e., two pairs of strongly degenerate heavy neutrinos. In this regime, they can 
only be discovered at LHC using pure Dirac neutrino signatures, which are more difficult 
to disentangle from the background. 

On the other hand, we find that in the complementary region of very small mixing 
e, e' <C 10~ 3 , m a R > 1 GeV, the lighter Majorana neutrinos ^4,^4 will lead to observable 
same-sign di- lepton signatures at LHC. A detailed study is missing, but previous results 
seem to indicate that a lower bound on 7714 of order 300 GeV could be set with 5 fb _1 of 
LHC data at 10 TeV [70]. 
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A. The neutrino mass two-loop integral 

The relevant integral is 

J = ,7(771,4, 771,4, mtt,"J,9,miy) = 
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ipq 

where 



V • Q 

1 I ( (A.l) 

((p + q) 2 - m|)((p + g) 2 - mf)(p 2 - m 2 )(q 2 - m|)(p 2 - m 2 v )(q 2 - mfy) ' 

d 4 p d g 



/"/ 



(2vr) 4 (2vr)^ 



We combine propagators with the same momentum by using 

1 1 [< dh 



then 



where 



(p 2 - m 2 a )(p 2 - m 2 w ) (m 2 - m 2 w ) (p 2 - t\ 

(pq) 

(p 2 -t 1 ) 2 (q 2 -t 2 ) 2 ((p + q) 2 -t 3 ) 2 ' 



J 



t Jpq 



ill r m l r m l r m l 

7 — 2 2~T7 — 2 TT7 — 2 2\ / ^ / dt2 / ^ 3 • 

Now we use the standard Feynman parametrization to combine the last two propagators 
and perform the integral in q, which leads to 

i i dx f p 2 



(4vr) 2 JJ l-xj p (p 2 - h) 2 (p 2 - t 3 /(l - x) - t 2 /x) 2 ■ 

The integral in p can be reduced by using an additional Feynman parameter and the final 
result can be written as 

(47r) 4 7 J y J t t 3 xy + t 2 (l-x)y + t 1 x(l-x)(l-y) ' 

The integrals in t\,t 2 ,t 3 can be done analytically and reduced to logarithms. The expres- 
sions obtained are complicated but can be used to feed the final numerical integration in 
x and y which converges smoothly for most of the parameters. The expression in (A. 2) is 
also very useful to obtain different approximations for small masses as compared with the 
largest mass in the integral. For that purpose one can use 

lim Jim -L_ / dt f(t) = lim f(a) = /(0) . 

a-tO b^a — a J a a^O 

Thus, for instance if 7714,777.4 ^> m a ,mp,mw ~ we can take t\,t 2 — > in the integrand 
and perform trivially the remaining integrals, 

1 1 777- 

J = ~ Ta VI 2 2 m 2' (A-3) 

(4tt) 4 m| - 7774 mi 

in agreement with the result in [48]. 

If 7713,7774,7774 3> m a ,mw ~ the integral can also be computed (take t\ — > in the 
integrand and perform the rest of the integrals). The result can be written in terms of the 
dilogarithm function Li2(x) and it is rather compact, 

r 2 m? 



J 



TT 



(4tt) 4 777 2 I 6 
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When mp — > it reduces, as it should, to ( |A.3| ) . 

We are especially interested in the case m a = mp = txie with rriE > mw, but 7774,7774 
could be larger or smaller than the (and even smaller than mw since we only know that 
7714 — m 4 > 63GeV). Some asymptotic expressions can be obtained when there are large 
hierarchies in masses 

1 1 , m\ 
(47r) 4 777^ 777y ' 

where mx is the heaviest of mj, mj, 7774, 777 and my the next to the heaviest of these 
masses. This expression can be used to perform analytical estimates but, since in the 
allowed range of masses the hierarchies cannot be huge, we do expect large corrections to 
these estimates. Fortunately, as commented above, the exact value of J for all values of the 
masses can be obtained numerically rather easily using ( |A.2| ). For fast estimates one can 
use 

1 1 , / 777? + 777 % \ 

J ~ I A \A ~2 2 2" m 2 ' 777 E , 7774, 777 4 > 777 W , 

(47TJ 4 ?7l| — 7774 — 777^ \ 77l| / 

which interpolates smoothly the different asymptotic expressions and reproduces the com- 
plete result with an error less than 50% in the worse case. 
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